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Abstract 

We study the Einstein field equations for spacetimes admitting a maximal 
two-dimensional abelian group of isometries acting orthogonally transitively on 
spacelike surfaces and, in addition, with at least one conformal Killing vector. 
The three-dimensional conformal group is restricted to the case when the two- 
dimensional abelian isometry subalgebra is an ideal and it is also assumed to 
act on non-null hypersurfaces (both, spacelike and timelike cases are studied). 
We consider both, diagonal and non-diagonal metrics and find all the perfect- 
fluid solutions under these assumptions (except those already known). We find 
four families of solutions, each one containing arbitrary parameters for which 
no differential equations remain to be integrated. We write the line-elements 
in a simplified form and perform a detailed study for each of these solutions, 
giving the kinematical quantities of the fluid velocity vector, the energy-density 
and pressure, values of the parameters for which the energy conditions are ful- 
filled everywhere, the Petrov type, the singularities in the spacetimes and the 
Friedmann-Lemaitre-Robertson- Walker metrics contained in each family. 
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1 Introduction 



One of the most successful ways of finding exact solutions of Einstein field equations 
has been the assumption of a certain degree of symmetry in the spacetime manifold 
[14fl . Although a few solutions are known without any symmetry at all, a systematic 



study of exact solutions with a given isometry group has only been performed when the 
dimension r of the isometry group G r is r > 2. In the particular case of G2 spacetimes, it 
is usually assumed that the isometry group is abelian and acts orthogonally transitively 
on 2-surfaces (that is to say, that the two planes orthogonal to the group orbits at each 
point are themselves surface- forming) . In the cosmological context, the two-dimensional 
orbits of the isometry group are taken spacelike (S2) and the matter content is usually 
assumed to be a perfect fluid. The Einstein field equations in this situation are still 
very complicated and further assumptions are needed to handle the problem. They 
may be very different in nature, for instance: separability of the metric coefficients, 
degenerate Petrov types for the Weyl tensor, particular equations of state, kinematical 
properties of the fluid velocity vector, Kerr-Schild ansatz, etc. Among these additional 
simplifications, we will assume the existence of conformal symmetries in the spacetime 
(for a definition of conformal symmetry see e.g. |14| and for a detailed study of some 
general properties of conformal Killing vectors see |l IJ). 

A conformal motion in the manifold is interesting for several reasons. First of all, it 
is a well-posed geometrical assumption which, therefore, does not depend on the coordi- 
nate system we are using. The Friedmann-Lemaitre-Robertson- Walker (FLRW) space- 
times (which are at present our best candidates to describe the real universe, despite 
them being highly symmetric) possess a fifteen-dimensional conformal group. Con- 
sequently, the inhomogeneous and anisotropic models containing a conformal Killing 
vector (CKV) are likely to contain FLRW limits and, therefore, are suitable as gen- 
eralizations of FLRW cosmologies (in particular, they can be used to test several per- 
turbative schemes used to study inhomogeneities in the universe). Another reason for 
studying inhomogeneous cosmologies admitting a CKV is the Ehlers-Geren-Sachs result 
[TOU which implies that observers measuring an isotropic microwave background radia- 



tion can exist only if the spacetime possesses a timelike CKV and the observer moves 
collinearly with it. 

Perfect-fluid spacetimes admitting a CKV have received considerable attention in 
the last few years. The case in which the CKV is inheriting (meaning that the con- 
formal motion maps the fluid world-lines into themselves) has been studied in ||. In 
particular, the spherically symmetric (|| and references therein) and the plane sym- 
metric spacetimes || admitting an inheriting CKV have been fully investigated. In 
it was shown that perfect-fluid spacetimes satisfying an equation of state p = p(p) and 
admitting a CKV parallel to the fluid velocity must be FLRW models. A first attempt 
to perform a systematic classification of spacetimes using its conformal isometry group 
has been performed in where some new exact solutions have also been found. All 
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this work showed that the perfect-fluid models admitting a CKV were rare but not 
impossible. In vacuum metrics, the most general spacetime admitting a CKV is 
either Minkowski or some type N solutions (representing certain plane wave solutions). 

The perfect-fluid case when the spacetime admits an abelian two-dimensional isom- 
etry group acting orthogonally transitively and one linearly independent proper con- 
formal motion (the homothetic case was studied in ||) has been considered only very 
recently |19], |§. In these articles, the authors classify the three-dimensional confor- 



mal group according to the possible inequivalent Lie algebras and find canonical line- 
elements for each case. Then, they concentrate on the diagonal metrics (which already 
excludes a number of Lie algebras, see below) and study systematically the possible 
solutions of Einstein's field equations for some selected Lie algebras (in our notation 
below, expression (pT3[), the Lie algebras considered in || correspond to Lie Algebra A 
with 6 = 0). 

In this paper we will complete this work by considering in full generality the case of 
spacetimes admitting an abelian orthogonally transitive G2 with spacelike orbits and 
one proper conformal Killing vector^ such that the three-dimensional conformal group 
has either timelike or spacelike orbits (for null orbits see ||19|| ). The final restriction we 
make is considering only the three-dimensional conformal Lie algebras for which the 
abelian isometry subalgebra is an ideal (i.e. the Lie bracket of the CKV with any of 
the two Killing vectors is a linear combination of only the two Killing vectors). We will 
not assume that the metric is diagonal (thus, we will study both classes B(i) and B(ii) 
in Wainwright's classification p0|). Obviously, we will assume b 7^ in expression (p~3|) 
below when restricting ourselves to the diagonal metrics. Although the assumption of 
two Killing vectors (KVs) and one proper CKV is the same as in [Tj|, |TJ]], in which 



the stationary and axisymmetric case was analyzed, the structure of the Lie algebras 
and differential equations to be solved turn out to be rather different. 

The plan of the paper is as follows. In section 2 we write down the Einstein field 
equations in non-comoving coordinates when the spacetime admits an abelian G2 acting 
orthogonally transitively on spacelike orbits. Then, we establish the necessary and suffi- 
cient conditions for a solution of the system of Einstein equations to represent a perfect 
fluid. We finish the section by showing that a new metric obtained by interchanging 
the coordinates t and x in all the metric functions of a given solution of the Einstein 
field equations has still the energy-momentum tensor with the structure of perfect fluid 
but with the 4- velocity switching from timelike to tachyonic and vice versa. This result 
allows us to simplify the work considerably as only one half of the cases need to be 
considered. In section 3 we write down the inequivalent three-dimensional conformal 
Lie algebras for which the isometry abelian subalgebra is an ideal. We then give the 
canonical line-elements for each inequivalent case (some of the forms of the Lie algebras 

1 Obviously, the conformal Killing vector is not unique since it is defined up to linear combinations 
of the Killing vectors and a non-vanishing constant factor. Throughout the paper we will refer to this 
equivalence class as 'the conformal Killing vector'. 
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and the line-elements are different to those in || although they are completely equiv- 
alent). In section 4 we study the non-diagonal metrics and proof that no perfect-fluid 
solutions with this specification are possible. Thus, in section 5 we concentrate on the 
diagonal metrics (when 6^0, see above) and find the general solution of Einstein's 
field equations. There exist four families of solutions (each of them containing arbi- 
trary parameters). These families are written explicitly in a simplified form in section 
6. Readers interested mainly in the resulting metrics may advance to this section where 
we give the Petrov type of the spacetimes, the form of the fluid velocity vector, the shear 
tensor components, the acceleration and the expansion of the fluid (vanishing rotation 
is already a consequence of imposing an orthogonally transitive Gz). The expressions 
for the energy-density and pressure of the fluid are written down and we find for which 
values of the parameters the energy conditions are fulfilled. We also discuss the exis- 
tence of an equation of state p = p(p), the singularity structure of the spacetime and 
the FLRW limit cases for each of the families. 



2 Some results on general abelian G<i orthogonally 
transitive perfect fluids 

As stated in the introduction, we are interested in cosmological abelian orthogonally 
transitive G2 perfect-fluid solutions of Einstein's field equations. It is well-known that 
for such spacetimes there exist coordinates {t, x, y, z} adapted to the Killing vectors in 
which the metric takes the form 



ds 2 ' 



dt 2 + dx 2 + F(t, x) (P _1 (t, x)dy 2 + P(t, x) (dz + W{t, x)dy)' 



S 2 (t,x) 

The two Killing vectors are obviously given by 

£ = d y , rf=d z . 

We are going to analyze the Einstein field equations for the metric (|l|) when the matter 
content corresponds to a perfect fluid. Thus, the energy- momentum tensor takes the 
standard form 



T a p = (p + p) u a Uf3 + pg a/3 , 

where p is the energy density of the fluid, p is the pressure and u a is the fluid 4- velocity. 
For orthogonally transitive G2 on S2 spacetimes, the fluid velocity vector is necessarily 
orthogonal at each point to the group orbits. Thus, the one-form u is always a linear 
combination of the coordinate one-forms dt and dx at each point and, consequently, 
the fluid flow is irrotational. In all the cases below, we will write the Einstein field 
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equations using orthonormal tetrads {0 a } such that 

b{t,x) b{t,x) 
Therefore, the fluid one-form will always take the form 

u = u 6° + u 1 1 , u 2 -ul = l, (2) 

so that the Einstein field equations, in units where c = 8ttG = 1, read 

S'oo = (p+p)ul -p, S 01 = (p + p) uqUu S n = (p + p)ul+p, 

S22 = S33 = p, S02 = S03 = S\2 = 5*13 = S23 = 0, 

where S a p stands for the components of the Einstein tensor in the {0 a } cobasis. 

Given that we are not assuming any particular equation of state for the fluid, these 
Einstein field equations can be rewritten in terms of the Einstein tensor only. The energy 
density, the pressure and the form of the fluid velocity vector can be found afterwards, 
once the field equations are solved. Due to the separation into two orthogonal blocks of 
the metric, some of the Einstein tensor components are identically vanishing. In fact, 
a direct calculation shows that S02, 5*03, 5*12 and 5*13 vanish and consequently, the only 
equations we need to consider are 

(S'oo + S22) (Sn — S22) — Sqi = 0, (3) 
S 22 - S 33 = 0, (4) 
S 23 = 0. (5) 

Equations (0)-(^|) are satisfied by all perfect-fluid metrics (with an orthogonally transi- 
tive G2) but not all the solutions of this system represent a perfect-fluid spacetime. Let 
us find which are the restrictions we need to impose to obtain the perfect-fluid solutions 
(in general these conditions can only be checked once the equations are solved). The 
quadratic equation (|3|) corresponds exactly to the vanishing of the determinant of the 
2x2 matrix 

I Sqo + S22 Sqi \ 
V S*oi Sn — S22 J ' 

which is clearly equivalent to the existence of two functions Uq and Ui such that 

Soo + S22 = gUq, Sqi = eUoUi, Su — S22 = tUf, 
where e is an appropriate sign. Only when the two functions Uq and Ui satisfy 

ul -U?>0 
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we can define the two components uq and u\ of the fluid velocity vector through the 
expressions 



ell = (p + p)u , eU 1 =(p + p)u 1 , with u - u t = 1, 

(which also define p + pasp + p = Soo — Su + 2S22). Obviously, for physically well- 
behaved perfect fluids e must be positive so that p + p > 0, but this is not necessary 
from a mathematical point of view for a solution of ©-© to be a perfect fluid. The 
separate expression for the density and pressure are obtained from S'22 = S33 = p. 
However, the solutions of (§)-(§) can also satisfy Uq — Uf < 0. When 



U — U 1 — Sqo + S22 — Sn — S : 



22 



the Einstein tensor (and therefore the energy-momentum tensor) takes the form 

S a/ 3 = K a Kp+ pg a0 

— * 

where K is a null vector lying in the two plane spanned by the vectors dt and d x . This 
kind of solution represents a radiative fluid with pressure (which should satisfy p = 
to be physically reasonable). As we are looking for perfect-fluid cosmologies we will not 
be interested in these radiative solutions here. 

Finally, when U$ - U 2 < 0, the energy-momentum corresponds to a perfect fluid 
but with the "fluid" velocity being spacelike. We will obviously discard this case in the 
following. 

To summarize, the conditions that the solutions of the system ©-(H) should satisfy 
in order to represent a perfect fluid can be written as 

S00 - S n + 2S 22 ± 0, — So ° Q +S "l Q > 0. (6) 

boo — On + ^022 

Our next aim will be to impose the existence of a proper CKV in the spacetime 
(one of the main assumptions in our work). Before entering into this we want to state 
an interesting property for general abelian orthogonally transitive G2 on S2 spacetimes 
which will simplify considerably our work. 



Lemma 1 If the spacetime 

ds2= SHt^) 

satisfies the equations (^), ^) and jj^), then the new metric 



dt 2 + dx 2 + F{t, x) (P~\t, x)dy 2 + P{t, x) [dz + W{t, x)dy) 2 )] (7) 



ds 2 



1 



S 2 (x,t) 



-dt 2 + dx 2 + F(x, t) (P _1 (x, t)dy 2 + P(x, t) (dz + W(x, t)dy) 2 ) 



obtained by interchanging x and t in all metric functions, also satisfies the set of dif- 
ferential equations (f|[) and ([|). 
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It is worth pointing out (and we will see that explicitly in the proof of the lemma) that 
if the first spacetime represents a perfect fluid then the second one is not a perfect-fluid 
solution because the conditions @ are not satisfied (the sign of u 2 is switched). Thus, 
this lemma does not allow us to find new perfect-fluid solutions from old ones. The 
interest of the result lies in the fact that the amount of work to be done is reduced 
considerably in many problems. In particular, we will see when imposing the existence 
of a CKV that we need to consider only one half of all different possibilities. Let us 
now prove the lemma. 

Proof: In order to avoid confusion, we will write the second metric (f|) using coor- 
dinates T and X instead of t and x. So the metric reads 



ds 2 



1 



S 2 (X, T) 



-dT 2 + dX 2 + F(X, T) (P~\X, T)dy 2 + 



P(X,T) {dz + W(X,T)dyf 



(9) 



Let us now consider the following line-element 



ds 2 



S 2 (t,x) 



adt 2 + adx 2 + F(t, x) P _1 (t, x)dy 2 + P(t, x) (dz + W(t, x)dy)' 



(10) 

where a — ±1. It is clear that when a = 1 this metric is exactly the same as (|7]) while 
when a = — 1 the coordinate change 



t = X, 



x 



T 



transforms the metric into (^). We are now going to analyze the Einstein tensor for the 
metric flT0|). We choose the tetrad 



S(t, x 



-dt, e 1 



S(t, x 



-dx, 



6> 2 



S(t,x) \ P(t 



F(t,x) 



,x 



dy, 6 s 



S(t,x) 



F(t, x)P(t, x)[dz + W{t, x)dy . (11 



The scalar products of these tetrad one-forms are, obviously, 



(d a ,O0) =diag(-a,a,l,l). 



The calculation of the Einstein tensor using this tetrad shows that in changing a — > —a 
the components S22, ^23, 5*33 switch sign whereas 5oo, £01, Sn remain unchanged (all the 
rest of components of the Einstein tensor are zero in both cases a = 1 and o = —1). We 
must now compare the components of the Einstein tensor of @ with the components 
of the Einstein tensor of (0) when in both cases an orthonormal tetrad with timelike 
6° is used (because this was used in (0) in order to write Einstein's equations in the 
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form ©-(§]))■ When a = 1 the tetrad fllT]) has already 6° timelike but when a = — 1 
we must interchange the superscripts to 1 because x is now the timelike coordinate. 
So, we find that the relationship between the Einstein tensors of the metrics (0) and 
(ID using respectively orthonormal tetrads with 6° being timelike is given by 

o(2) _ q(l) q(2) _ Q (l) q(2) _ q(l) 

°00 — D ll ) °01 — °01 ; D ll — °00 5 

q(2) _ q (l) q (2) _ q (l) q (2) _ Ml) 

°22 — °22 ) °23 — °23 > °33 — °33 j 

where the superscript (1) denotes the first metric ([?[) and the superscript (2) denotes 
the second metric Consequently, if S^p satisfies the equations ([|), (|J) and (|5|) 

(2) 

then S a p also satisfies the same equations and the lemma is proven. Recalling that the 
conditions for the solutions of these equations to represent perfect fluids are (H) and 
given that 

q (2) q (2) _ q(l) q (l) q (2) q (2) _ q (l) q (l) 

°00 i °22 — °11 J 22 i °11 °22 _ °00 i" °22 ) 

we find that in spacetime regions where one metric represents a perfect fluid, the other 
represents a tachyon fluid and vice versa. 

□ 

Until now the discussion has dealt with general abelian G2 orthogonally transitive 
perfect fluids. In the next section we will impose the existence of a CKV, classify the 
inequivalent Lie algebras and write down canonical line-elements for each resulting case. 



3 Inequivalent Lie algebras and line- elements 



The main assumption in this paper is the existence of a proper conformal motion in the 
manifold. We do not assume any further conformal symmetries in the manifold, but we 
assume that the two abelian KVs and the CKV, k, span a three-dimensional conformal 
Lie algebra which, in general, satisfies 



0. 



6, k = 7i£ + 7 2 r/ + 73/c, ff, k = 5i£ + 5 2 f] + 8 3 k 



(12) 



where all 7's and <5's are obviously constants (in the presence of a fourth symmetry 
the existence of the three-dimensional Lie subalgebra (0) could be violated). As a 
short calculation shows, all Lie algebras (|T2"D are solvable (see e.g. |18j for a definition 
of solvability). We in the following make a stronger assumption. We assume that 
the two-dimensional abelian subalgebra generated by the two KVs in ( |T2"D is an ideal 
(understanding the multiplication in the algebra as taking the commutator), i.e. we 
assume 73 = £3 = 0. The line-element and consequently the Einstein field equations 
take very different forms whether the isometry subalgebra is an ideal or not. Thus, 
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the two cases must be treated separately. In this paper we will concentrate only on 
the case when the isometry subalgebra is an ideal. Thus, we assume from now on 
that 73 = #3 = 0. We must now find canonical forms for the inequivalent Lie algebras 
contained in fll2j) for 73 = £3 = 0. In order to do so, we can exploit the freedom in 
performing linear transformations of the KVs £ and ff. Under these transformations, 
the 2x2 matrix 

( 7i 72 \ 
\S 1 5 2 J 

transforms like an endomorphism. By doing this, we can reach a canonical Jordan form 
for this matrix. Three different possibilities arise depending on whether the matrix 
diagonalizes with real eigenvalues (Case A), diagonalizes with complex eigenvalues (Case 
VII) or does not diagonalize (Case B). We used the term VII because that Lie algebra 
corresponds to the Bianchi type VII in Bianchi's classification of three-dimensional Lie 
algebras. Each one of the other two cases, A and B, contains several different Bianchi 
types and no similar Bianchi type name can be used (I, III, V, VI are in A and II, IV are 
in B). In the complex case, we then perform a transformation to obtain a real valued 
matrix. The three possibilities read, finally, 



Lie Algebra A 



v] = 0, [£, k\ = - (c + b) £ [ff, k\ = - (c - b) ff, 
where b and c are arbitrary (possibly vanishing) constants. 

Lie Algebra B 



(13) 



0. 



£> k \ = + ar ^' [ft k 



where a is a no n- vanishing constant. 
Lie Algebra VII 



e,^J=o, [tk 

where a and c are constants (a ^ 0). 



\ c £- a V> [if, k] = at[+ -cff. 



In order to find the line-elements for each of these Lie algebras we must solve the 
conformal Killing equations 
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where represents the Lie derivative along the vector field k and $ is a scalar function 
usually called scale factor. Given that we will concentrate on spacetimes admitting a 
proper CKV, the scale factor will be assumed to be non-constant (and in particular not 
identically zero). When solving the conformal Killing equations, it is found that three 
different cases must be considered depending on the character (spacelike, timelike or 
null) of the orbits of the three-dimensional conformal group (which is generated by the 
Lie algebra {£,f),k}). In this paper we will not study the case when these orbits are 
null and, instead, we will concentrate on the cases when the orbits are either timelike 
or spacelike. Let us first write down explicitly the canonical line-elements for each of 
the Lie algebras A,B and VII when the orbits of the conformal group are timelike. 

Lie Algebra A 



ds 2 



S 2 (t,x) 



-dt 2 + dx 2 + F(x)p-\x)e- { - h+c)t dy 2 + 
+ F(x)P(x)e (b ~ c)t (dz + W(x)e- bt dy) 



The expression for the conformal Killing vector in this metric is 
The orthonormal tetrad we will use for this Lie algebra is 



(14) 



6° = \dt, 1 = ^dx, d^^Me-^dy, 6 s = ^-e^ (dz + We' bt dy) .(15) 



Lie Algebra B 



ds 2 



S 2 (t,x) 



dt 2 + dx 2 + 



+ F(x)e- Ct (P- I (x)dy 2 + P{x) (dz + [W(x) + at] dy) 2 ) 



(16) 



The conformal Killing vector in this metric reads 

k = d t + ^cyd y + (ay + ^cz^j d z . 
The orthonormal tetrad we will use is 

0° = -dt, 1 = -dx, 



e 2 = 



1 F 

sVp ( 



£ 2 f dy, s = 



VFP 



e~i* (dz + [W + at] dy) . 



(17) 



10 



Lie Algebra VII 



In this case the explicit form of the metric is much longer than before. We will give, 
instead, the explicit form of the orthonormal tetrad we will use in our calculations. It 
reads 

00 = I rft, 1 = - dx, 

S(t, x) ' S(t, x) 



6> : 



S(t,x)\ P(x) 



Fix) 



e 2 * [cos (at) dy — sin (at) dz] , 



3 = \ F ^ P ^ e i [ cog ^ _ w ^ gin ( at ^ dz + [ gin ( at ) + w ^ cog ^ dy \ _ 



S(x,t) 

The conformal Killing vector in this metric is 



k = d t + ( ^cy + az^j d y + (^-ay + ^cz ) d z . 



The scale factor of the CKV for each one of these metrics reads 

*=-§?■ (19) 

For the remaining cases in which the orbits of the conformal group are spacelike, the 
line-elements can be obtained directly from the previous ones by simply interchanging t 
and x in the metric coefficients as described in section |[ It is now obvious that making 
use of the lemma 1 we can restrict the study to the case when the orbits of the conformal 
group are timelike because in finding the general solution of the three equations @, (f|) 
and (|5|) we will also find the general solution for the metrics obtained by interchanging 
the roles of x and t in the metric coefficients, which correspond exactly to the metrics 
when the orbits of the three-dimensional conformal group are spacelike. 

In the following section we will concentrate on the non-diagonal metrics and will 
prove that there do not exist perfect-fluid solutions under the assumptions of this paper 
when the line-element is non-diagonal. 



4 Non-diagonal metrics 

In this section we are going to prove that all the perfect-fluid solutions of Einstein field 
equations for the spacetimes we are considering in this paper (see beginning of section 
3 for a clear account of our assumptions) must be diagonal. In the previous section we 
found the canonical line-elements for the three possible Lie algebras (A, B and VII). 
It is clear from the expressions for these line-elements that all the metrics belonging 
to Lie algebras B and VII must be necessarily non-diagonal (irrespective of any field 
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equations) and that only for Lie algebra A the diagonal case is possible (setting W = 
in flliD ). Thus, in order to prove that no non-diagonal perfect fluids exist, we must 
show that the metrics in Lie algebras B and VII do not admit perfect-fluid solutions at 
all and that only the diagonal metrics in Lie algebra A admit perfect-fluid solutions. 
These results will be stated in the following three propositions (each one corresponding 
to one of the Lie algebras A,B and VII). Let us start with the simplest case, namely 
the one corresponding to Lie algebra B. 

Proposition 1 No perfect-fluid solutions exist for the conformal Lie Algebra B with k 
being a proper CKV. 

Proof: We will concentrate on the two Einstein field equations (f|) and (||) which, 
using the tetrad ([H]), read respectively 

^-W' 2 P 2 + a 2 P 2 - 2^ = (20) 
,<7 \ 1 1 ,<7 , 

= (21) 

where the prime denotes derivative with respect to x and the comma indicates partial 
derivative. Let us first consider the situation in which P' ^ 0. The first equation above 
shows that the function S must take the form of a product of a function of x and a 
function of t, S = G(x)H(t). Then, the second equation reads 

which immediately implies that H/H must be a constant (remember that for the Lie 
algebra B the constant a cannot vanish) and then the expression for the scale factor 
(|19|) implies that the CKV is in fact homothetic. Thus no perfect fluid with a proper 
conformal motion does exist in this case. It only remains to analyze the situation when 
P satisfies P 1 = 0. The equation (p0|) gives simply 

W = eia, 

where t\ = ±1. The second equation ( |2"TD takes the form 

IF' 1 S t S x 

which can be integrated to give 

S = y/Fexp ( e l^ x J U(t+ e-yx) 
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with U an arbitrary function of t + ex. The quadratic Einstein equation ([3]) takes now 
the very simple form 

/ jpii 771/2 \ 2 

1 = ^ F = F^, 



2 



with constants Fq and (3. Now we can evaluate the components of the Einstein tensor 
to find 

•SoO = 5*11 = CxSoi, 5*22 — 5*33 = S*23 = 0, 

which shows that the energy-momentum content is that of a null fluid without pressure. 
Thus, we can conclude that no perfect-fluid solutions with a proper CKV exist for the 
Lie algebra of the type B. Let us emphasize once more that this analysis covers both 
cases, when the orbits of the conformal group are timelike and spacelike. 

Let us now prove a similar result for Lie algebra VII, i.e. that no perfect-fluid 
solutions arise for Lie algebra VII. 

Proposition 2 No perfect- fluid solutions of Einstein's field equations do exist for space- 
times 

- admitting a three-dimensional Lie algebra of class VII of two KVs and one proper 
CKV with 

- the two-dimensional isometry group being maximal. 

Proof: The proof for this case is similar to the previous one and we will only outline 
the main steps. We concentrate, as before, on the two Einstein equations (|) and (||) 
which, using the tetrad fll8|) , read respectively 

F'P' ( P'\ P'S X Tir S t 

h — - 2 - AaW— - 

F P \ P j PS S 



a 2 



-W' 2 P 2 + a 2 P 2 (W 2 + l) --^-2acl^ = (22) 
, / 1 F' P' S x \ (S t c\ f Tjr2 , 1 

w {2F + p-1r ) +a {-j + 2) r 2 + 1 -p2- 

+1.W" + a 2 W (w 2 + 1 + j^j = 0. (23) 

We must again distinguish between two cases depending on whether P' = or not. 
When P = const, the only non-homothetic solution of the system (f2~2j)-(|23|) is P = 1 
and W — 0. Then, the tetrad reduces to 

dt dx 1 

0° = —, 1 = -z-, 3 = -VFe-^ (cos {at) dy -sin {at) dz ), 
o o o 

4 = ^VFe'^ (sin (at) dy + cos (at) dz ) 
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and the line-element is simply 
ds 2 ' 



s(x,ty 



-dt 2 + dx 2 + F(x)e~ ct (dy 2 + dz 



which has a three-dimensional isometry group (against our assumptions). Thus, we 
must consider the case P' 7^ 0. It is not difficult to see that the only possibility for 
proper conformal solutions is 

S = G(x)H(t + aK(x)), (with K' = -2Wy)j 
P' G' 1 (F' \ , m 2 1 

where m is a constant, e 2 = ±1 and H is an arbitrary function of its argument. Imposing 
now the quadratic Einstein equation (|3]) we simply find 

where Fo and (5 are constants of integration. As in the previous case, the evaluation of 
the Einstein tensor for this solution gives us 

5*00 — <Sll — — ^'S'oij 5*22 — ^33 = S23 = 

so that the matter content is a null fluid without pressure and the proof of the propo- 
sition is completed. 

Only Lie algebra A remains to be analyzed. Let us prove that all the perfect-fluid 
solutions with a proper CKV belonging to Lie algebra A must be diagonal. 

Proposition 3 No perfect-fluid solutions of Einsteins field equations do exist for space- 
times 

- admitting a three-dimensional Lie algebra of class A of two Killing vectors and 
one proper CKV with 

- the maximal two-dimensional isometry group being maximal and 

- the metric being non- diagonal. 

The proof is again similar and we will just sketch it. The two equations (f|) and (||) 
read respectively 

+(?)'" 2 T S 1 + 26 T " + + bc = °> (24) 

W>^ + y-^)+\w" + h(b-c)W-bW^ = 0. (25) 
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The P' = subcase gives rise either to homothetic spacetimes or to metrics with 
a three-dimensional group of isometries acting on two-dimensional spacelike surfaces. 
Thus, we can move to the general case P' 7^ 0. Equation (p^ ) gives 

S = G(x)H(t + bK{x)) with K' = y r 
which inserted into (|25|) produces 

w ' (It + t - %) + \ w " +1 2 b{b - c) w = I" ( w + w t) ■ 

The conditions that the CKV is proper implies that both sides in this equation must 
vanish. In particular, the right-hand side implies that either b = or W = a/P (where 
a is an arbitrary constant). Let us start with 6 = 0. The two equations (|24]) and (|25|) 
take now the form 

and the second one can be integrated to give 

C 2 

W = Wn—— 

where Wo is an arbitrary constant. The case Wq = gives a diagonal metric and, 
therefore, we can assume Wo ^ 0. Then, equation (^) can be integrated to give 

±- 2 =6 2 -(W + [3f (28) 

(where /3 and S are arbitrary constants of integration). The metric is apparently non- 
diagonal, but a trivial calculation using (|28p shows that the linear coordinate change 
in the block {y, z} 

y = Y + Z, z = {(3 + 5)Y + {(3-5)Z 

brings the metric into a diagonal form and therefore no proper non-diagonal solutions 
do exist in this subcase (notice that this coordinate change is singular only when 5 = 
which is impossible due to (|28|)). 

We can now consider that last possibility, namely, when b 7^ and W = a/P (with 
a 7^ in order to consider non-diagonal metrics). Substituting W everywhere, it is easy 



(26) 
(27) 
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to see that the general solution of the full set of Einstein field equations ([!]), (f|) and 
©is 

P' G' 1 . , F' 

-p =(Tb > -q = 2 ( n + ac > ' y = n ' 

where a = ±1 and n is an arbitrary constant. As in the two previous cases this solution 
corresponds to a null fluid and the proof of the proposition is completed. 

To summarize, we have proven in this section that there do not exist any non- 
diagonal perfect-fluid solutions for spacetimes with 

- a maximal abelian orthogonally transitive G<i on Si 

- admitting one proper CKV (with conformal orbits either spacelike or timelike) 
such that 

- the subalgebra generated by the two Killing vectors is an ideal in the Lie algebra 
spanned by the two KVs and the CKV. 

In the following section we will concentrate on the diagonal case and find all the perfect- 
fluid solutions (with b ^ 0, see the introduction). 



5 Diagonal perfect-fluid solutions 

We know that only the Lie algebra A admits diagonal metrics (all metrics in Lie algebras 
B and VII are necessarily non-diagonal and, as we have seen in the previous section, 
none of them satisfies the perfect-fluid Einstein field equations). In order to restrict the 
metrics for Lie algebra A to be diagonal we must set W = in (|i~4f) (and also in the 
corresponding tetrad (|T5D). Now, the calculation of the Einstein tensor in this tetrad 
gives S23 = so that the Einstein field equation (|]) is identically satisfied (this can be 
seen directly from after substituting W = 0). Equation fl2~3| ) with W = gives 
(P = const, would imply 6 = and a third KV): 

S = G(x)H(t + bK(x)) with K' = —, 

be = 0. (29) 

Now, we must inevitably analyze the more difficult quadratic equation @ which takes 
the form 

E x (x) + |s 2 (x) + ^S 3 (x) + f§S 4 (x) + |s 5 (x) = 0, (30) 







-4- 
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where Sj(x) (i — 1 to 5) are expressions containing P, F and G and their derivatives 
and which depend only on the variable x. The dot means from now on derivative 
with respect to the variable t + bK(x). In the Appendix we find all the different 
possibilities compatible with this mixed equation containing both independent variables 
x and t + bK. Let us here just summarize the results. 

Case 1 



£ x = £ 2 = £ 3 = £ 4 = £ 5 = 0, H arbitrary 



Case 2 



£4 = 0, £1 = WI1E5, £ 2 = m 2 £ 5 , £ 3 =m 3 £ 5 , E 5 ^0 

H H 2 H 

™\ + m 2 — + m 3 — + = 0, 

where mi, m 2 and m 3 are arbitrary constants. 
Case 3 



£1 = fciE 4 , £ 2 = fc 2 £ 4 , £ 3 = £; 3 £ 4 , £ 5 = /c 5 £ 4 , S 4 ^0 
, H , H H , H 

where k±, k 2 , k 3 and k 5 are arbitrary constants. 
Case 4 

£4 

£1 = kiE 5 , £ 2 = fci£ 4 + fc 3 £ 5 , £ 3 = /c 3 £ 4 , £5 7^ 0, — ^ const. 



£ 



5 



H H 

kl + h H + H = ° 

where k\ and fc 3 are arbitrary constants. 

The explicit expressions for £j will be necessary later to study these four different 
cases. They can be directly obtained from the following expressions. 

£1 = Z X V X - Wl £ 2 = Z X V 2 + Z 2 V 1 - 2W 1 W 2 , £ 3 = Z 2 V 2 - W*, 
£ 4 = -2 (2W 2 R + Z 2 R 2 + V 2 ) , £ 5 = -2 (2W 1 R + Z X R 2 + Vi) , 
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where Z\, Z%, Vi, V2, W±, W2 and R are given explicitly in terms of the metric coefficients 
by 



1 2 F 



Z 1 = 


IF" 
2~F ~ 


f' a 
f"g 


+^ 


Z 2 EE 


F' P 

— c — , 

P P'' 


IF' 2 
2 F 2 " + 


f' a 

f"g ~ 


2^- 
G 


ip' 2 

f 2 pa + 


l -c 2 
2 ' 


P 2 

V a = c- 2b 2 c— 2 


Wi = 


cF' 

Yf ~ 


bP' 

2 p~' 


^2 E 


= 2^ 
G' 


P 



p> P 
} ~F~p~r 



As explained in the introduction, the case when the constant b in the Lie algebra A is 
identically zero has been exhaustively analyzed in 0. Thus, we can restrict our study 
to the case when b 7^ 0. We have to solve the four different sets of differential equations 
Case 1 to Case 4. Obviously, the two first cases are simpler because of the equation 
£4 = 0. Let us first concentrate on these two cases. It is convenient to define three new 
functions q(x), l{x) and m(x) through the relations 



P' , F' 



a 



i)q, 7? = - c + t- 



G q 



in 



The equation ( p9|) reads now 



q' +(c + I) q 2 



m 
~2 







which allows us to obtain m in terms of q and I. Using this expression for m, the 
equation S 4 = takes the form 



2q' + + 2c) (l - q 2 ) = 0. 



(31) 



The solution q = ±1 of this equation also implies S 5 = and, therefore, the only 
possibility is Case 1. The expressions for S 2 and S 3 are identically zero while the 
vanishing of imposes I' = 0. Now all the equations are satisfied. However, evaluating 
the Einstein tensor for this solution we find 



5*00 + S; 
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S11 — S; 



22 



±S t 



01) 



S- 
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S: 



.33 







so that it does not represent a perfect fluid. Thus, we can assume q 2 7^ 1 and obtain I 
from (RID as 



I 



-2c 



2q' 



q 2 - 1 

The expressions for S 5 and S 3 read now, respectively, 



S n = — 



1 



<? 4 

S 3 



cq 



1 ( 1 1 2 

- 2 (q + cq - 
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Let us start by analyzing Case 1. Given that both £5 and E3 must vanish we immedi- 
ately find 



q' = c (l - q 2 

where e is a sign. Now, all £j are identically zero and we have another solution of the 
Einstein field equations. The calculation of the Einstein tensor shows that the matter 
content is indeed a perfect fluid. However, the spacetime is conformally flat and both 
the density and pressure depend only on the variable t + bK. Thus they satisfy an 
equation of state p = p(p) and, consequently, the spacetime is a Friedman-Lemaitre- 
Robertson- Walker cosmology. Consequently, Case 1 does not give new perfect-fluid 
solutions and we can move on to Case 2. 
The equation £3 — m 3 £ 5 = reads 

(1 - m 3 ) (q + cq 2 - c) 2 = m 3 (b 2 - c 2 ) (q 2 - if . (32) 
When m 3 7^ 1 this equation immediately implies 

q > = n (l - q 2 ) , (33) 



where n is a constant. When m 3 = 1 the equation ( p2|) implies b = ec and then, 
combining the other two differential equations Si — miEs = and £2 — = (two 
equations for one unknown q(x)), it is not difficult to show that the relation fl3~3|) must 
still hold. Substituting this expression for q' into all Ej we find that all the equations 
of Case 2 are fulfilled, and they simply give m 1; m 2 and m 3 in terms of the constants 
b, c and n. The explicit expressions are 

(c — n) 2 (n — c) 3 + 3c (n — c) 2 + c (b 2 — c 2 ) 



m 3 = — — , 7B2 = 

n 2 — 2cn + b 2 n 2 — 2cn + b 2 

1 (c 2 - n 2 ) (c 2 - 2cn + b 2 ) 
mi = ~- 



a 



2 



2cn + b 2 



and the metric functions are given by 



— = bq{x), — = (2n - c) q(x), ~q = ~ c > ^ + 2~q(x) ) Q =n V~ ( l 

where b, c and n are arbitrary constants. Let us label this solution as Solution A. 
We will write down the explicit line-element and analyze the physical content of this 
solution in the following section. 

Thus, we have exhausted the Cases 1 and 2 and we have to deal now with the much 
more complicated Cases 3 and 4. We have to take advantage that in both cases the 
equation S 3 — /C3S4 = holds. Therefore, we will try to write this equation in the 
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simplest possible form so that the problem becomes more tractable. To that aim let us 
define three functions r(x), u(x) and v(x) through the relations 

F' , , , , G' , s P' 



— = 2r(x) +u(x), — = r(x), — = Kb^/v(x), 

where k is the sign of P' / (bP) (so that the square root is non-negative definite, as 
usual). Now, equation fl29|) allows us to obtain u in terms of v as 

2ckJv + v' 

u = 

2v 

and substitute it in all other expressions so that only the two functions v(x) and r(x) 
remain in the equations. 

Let us consider the situation when v is a constant. Then we must assume o / 1 
because otherwise we would have £4 = and would be in a previous case. Now, the 
equation S3 — ^3X4 = immediately implies r = const. So, we have 

f' c a p' 

— = 2ai , — = ai , — = bao 

F a 2 G P 

(a\ and ct2 constants) and we certainly have a solution of Case 3 (Case 4 is impossible in 
this situation because E5 is constant and therefore proportional to £4). The equations 
simply give k\, &2, k 3 and k§ in terms of the constants c, b, a\ and a-i- This is a solution 
of the Einstein field equations and its matter content is a perfect fluid, but it is trivial 
to checked that the metric contains a homothetic Killing vector (besides the proper 
CKV) and therefore does not belong to the class we are interested in this paper. 

So, we can assume from now on that v is non-constant and perform the following 
change of variables 

dv = 2k\/v (v — 1) \x( v ) ~ ^3] dx 

so that now the independent variable is v and x( v ) is the unknown. Let us also define 
a function d(v) through 

Ik 

r = -—= (u - x + 2k 3 ) 
2 \/v 



so that we replace r by u everywhere. With all these changes and redefinitions, the 
common equation E 3 — /C3S4 = takes the simple form 



X*-(c-k 3 y\ =u 2 -(c-k 3 ) 2 . (34) 

This equation is linear in the independent variable v and, therefore, it would be conve- 
nient to consider x as the independent variable and v(x) the unknown. This change of 
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variables can be done unless x is constant. Thus, we must consider the two cases x — 
const, and x non constant separately. Let us first assume x — const. 

Now, we have only one unknown u(v ) which must satisfy an overdetermined system 
of differential equations (four for Case 3 and three for Case 4). It is not difficult to find 
the solutions which are: 

Case 3 and x = const. 

The solution in this case is 

X = k 3 - aib, u = k 3 - (Txb, c = (nb, 

3b 2 b 3 

k x = — — (<7i&-A&), h = - {-b + Aaxka) , h = -a x b, 

where a x — ±1 and b and k 3 are arbitrary constants. Let us label this solution as 
Solution B for further study in the next section. 



Case 4 and x = const. 



Now, two different solutions are possible. The first one is 

X = c — k 3 , u = c — k 3 , k\ = 0, b 2 + c 2 — 2ck 3 = 0, b, c arbitrary 

which again will be studied in the following section. Let us call it Solution C. The 
other possibility reads 

c 2 -b 2 

uo = a 2 (c - aib) y/v, p = c - crib, k 3 = c, k x = , 



where a 2 is also a sign. This solution, however, can be seen to represent a null fluid 
(with non-zero pressure, in general) and, therefore, will not be analyzed any further. 

Having written all the solutions of the equations when x = const, let us consider 
the more general case when x is non-constant. Now, as stated above, it is convenient 
to perform another change of variables and consider x as the independent variable and 
v and ijj as unknown functions depending on x- Equation S 3 — k 3 T>i = 0, (|34"D, allows 
us to obtain v in terms of x as 

u 2 -(c-k 3 ) 2 

V = 2 

X 2 - (c - k 3 ) 

and to substitute this expression in all the Sj. Now we have only one unknown u)(x) 
which must satisfy an overdetermined system of differential equations (three differential 
equations for Case 3 and two for Case 4). The investigation of the compatibility of 
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these systems of differential equations is now rather lengthy and makes essential use 
of computer algebra | 13| . Special care is needed to combine the equations in the right 



order to avoid too long expressions. Let us here simply summarize which are the results 
obtained. It turns out that there exists one family of solutions for Case 3 and another 
family for Case 4. They are, explicitly, 

Case 3 and x non constant 

uj = x — 2o"i& + 2k 3 , c = <Tib, ki = 0, 
/c 2 = - (2cri&;3 — b) , k 5 = k 3 b and k ?J arbitrary. 

This solution will be labeled as Solution D and will be studied in the following section. 

Case 4 and x non constant 



W = X 1 + 



k\ - 4*4 - b 2 



\ (x-h + b)( x -h-b) 



k\ and k 3 arbitrary. 



The Einstein tensor for this solution satisfies S00+S22 = ^ii — S22 = — £01 and, therefore, 
as discussed above, represents a null fluid. Consequently, it will not be analyzed any 
further. 

To summarize, we have found four families of solutions which represent perfect fluids 
when the spacetime is assumed to have an abelian orthogonally transitive G2 on S2 and 
also one proper CKV which spans together with the two KVs a Lie algebra generating 
either timelike or spacelike orbits such that the isometry subalgebra is an ideal. These 
families have been labeled as Solution A, B, C and D respectively and will be studied 
in the next section. 



6 Analysis of the Perfect-Fluid Solutions 

The aim of this section is to write down the line-elements for the four solutions found 
in the previous section (Solutions A, B, C and D). In order to do so, we must integrate 
back some changes of variables which helped us to solve the Einstein field equations. 
We will not give the details as the calculations are not difficult and can be easily 
reproduced. In some cases we have also performed some redefinitions of constants and 
some coordinate changes in the block {t, x} as well as rescalings of the coordinates y 
and z in order to write the metrics in a simplified form. Therefore constant names used 
in this section are unrelated to previous appearances in earlier sections. After giving 
explicitly the line-element (no differential equations remain to be solved) we will study 
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some of their properties. In particular, we will give the Petrov type, the fluid 4-velocity, 
the kinematical quantities, the energy- density and the pressure. 

We will also study in which spacetime region the solution represents a perfect fluid. 
As discussed earlier, there is a symmetry between exchanging x and t in the metric 
coefficients and switching between a timelike perfect fluid and a tachyonic fluid. In the 
case of metrics A and D there are regions of both, perfect and tachyonic fluid. Here a 
switch x <-> t might be of interest to exchange the physically meaningful regions. 

Further, we describe for which ranges of the parameters the spacetime satisfies 
energy conditions everywhere (in the perfect-fluid region) and whether the fluid satisfies 
an equation of state p = p(p) or not. 

A brief description of the curvature singularities of the spacetime and an identifi- 
cation of the spatially homogeneous and isotropic limits in each family are also given. 
All four families have the magnetic part of the Weyl tensor with respect to the fluid 
4-velocity being non-zero. 

Tensor components have been calculated with the help of the computer program 
CLASSI ]T5| using the orthonormal tetrad 



6 a = y | g aa \dx a , (no summation over a) 

for each metric (all four line-elements are diagonal). 



6.1 Analysis of the solution A 

After rescaling the coordinates {t,x,y,z} and a coordinate change e*cosh(x) — > t, 
e* sinh(x) — > x the metric can be written as 

ds 2 = S- 2 (-dt 2 + dx 2 + t l - a ~ b dy 2 + t l - a+b dz 2 ) 

where S = t- a / 2 x 1 ~ a2 / q and a, b are essential constants that are free apart from the 
condition, that q defined as q = b 2 + a 2 — 2a — 1 has to be non-zero. Coordinates t, x 
are restricted to < t and < x to allow arbitrary real exponents in S. We will see 
below that in order to satisfy positivity conditions of energy, the parameter a has to 
take values for which there is a curvature singularity at t = such that the restriction 
t > is no loss of generality. As a switch b <-> — b is equivalent to y <-> z, we consider 
only b > 0. The metric has two conformal Killing vectors 

ki = d x , k 2 = 2td t + 2xd x + (1 + a + b)yd y + (1 + a — b)zd z . 

With the further abbreviation D = q 2 x 2 — Aa 2 t 2 the fluid velocity one-form has com- 
ponents u a = — sign(g)D _1//2 (gx, 2at, 0, 0). As announced before, lower Greek indices 
denote components in the frame 6 a . An overall minus sign was chosen in order to 
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have a future pointing u. To exclude null or tachyonic fluid the coordinates have to be 
restricted to guarantee D > 0, i.e. 



- > 2 
t 



(35) 



The Petrov type of the metric is / and the non-vanishing Weyl spinor components \l/ 
are 

abS 2 _ (6 2 + a - l)S 2 



I2t 2 



In the case b = 0, a = 1 the metric is conformally flat. Then the metric has 6 KVs, 
shows p — p < and obeys an equation of state which is 



P + P 



8 



P-P 
12 



- 1 



It therefore is a special FLRW spacetime. Otherwise, in any one of the four cases 

a = 0, 6 = 0, 6=1, 6= | 1 - a \ 

the metric is of Petrov type D. The expansion of the fluid is 

S 



e 



[48(a 2 - gKt 4 - lQa(a 2 - q)q 2 t 2 x 2 + (a + 2)gV]. 



2\q\txD 3 / 2 

The acceleration and the acceleration scalar of the fluid are 

(2a + q)q 2 x 2 S , 
a a = — 2 (2at, qx, 0, 0) , 

a a a a = (2a + q) 2 q 4 x 4 S 2 /D 3 . 
The non-vanishing components of the shear tensor a a p read 



o"oo 



2at 
qx 



2+2 



coi 



4aH 
q 2 x 2 



an = - 



Aa 2 \q\txS[D(l -a)- Aat 2 (2a + q)\ 



0"22 



0"33 = 



3D 5 / 2 

\q\xS[(l -a- 3b)D - 4(2a + q)at 2 ] 
6tDV 2 

\q\xS[(l - a + 3b)D - 4(2a + q)at 2 ] 
6tD 3 / 2 



The shear scalar is 



q 2 x 2 S 2 



(3b 2 D 2 + [Aat 2 {a 2 + 6 2 - 1) + (a - l)D} 2 ) . 
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With the energy density 



and the pressure 



P 



(a 2 -q)S 2 [q 2 x 2 - 12(a 2 - q)t 2 ] 
Aq 2 t 2 x 2 

(a 2 - q)S 2 [q 2 x 2 + A(a 2 - 3q)t 2 ] 
Aq 2 t 2 x 2 



we obtain 



(a 2 -q)S 2 D 2(a 2 -q)(2>q-2a 2 )S 2 

P + P= — n 9.9 9 — ' P-P 



2q 2 t 2 x 2 ' q 2 x 2 

To have the last two quantities non-negative and therefore p non-negative, a and b have 
to satisfy 

a>-l/2, 2a+l>b 2 >2a + l-a 2 /3. 

For a > —2 there is a big-bang like singularity at a finite distance. There are two 
possibilities to have an equation of state p = p(p): 

- For q = —2a both, p and p, are functions of x/t. 

- For q = 2a 2 /3 we have p = p= (a 2 x 2 - 9t 2 )/(12t a+2 x 3 ). 

The metric A has an abelian three dimensional conformal subgroup and therefore 
belongs to a class of metrics discussed previously by A. M. Sintes in and A. M. Sintes, 
J. Carot and A. A. Coley in |3] although no explicit expressions for the metric were given 
there. In Q it is explained that the energy conditions p + p > 0, p — p > can only be 
fulfilled in a special subclass to which metric A does not belong to. In contradiction we 
found that for any a > — 1/2 there exists a range of values for b such that both energy 
conditions are satisfied in the complete perfect-fluid range. 



6.2 Analysis of the Solution B 

The metric is 



dt 2 



sinh 



2+2m, 



+ 



dx 2 



x 



sinh 



2+2m, 



+ 



X) 



e -2t 

-dy 2 + cotanh 2 (a;)(i,2 2 



sinh (x) 



where S depends on the variable u = t + log sinh (x). This coordinate approaches 
asymptotically a null coordinate for big x. The conformal Killing vector of this metric 



is 



k = d t + yd y . 
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We define a new function z{u) as 



dS 
du 



— = (z-l)S. 



The ordinary differential equation that S(u) had to satisfy has in general two constants 
of integration. However, by adding appropriate constants to the coordinates t and 
x and rescaling all the coordinates, it can be easily seen that there only appear two 
inequivalent cases (we are discarding the homothetic solutions because they have been 
already investigated in H). The explicit expressions for S and z are: 



l+m. 



S = be-— u cosh( 



-u 



1 — m + vl + m 2 tanh 



Vl + m 2 



-u 



S = be 2 M s inh( 



«) 

1 — m + Vl + m 2 cotanh 



where b is an arbitrary constant. So, the metric has two essential parameters: the 
constant b gives a global scale to the metric and the constant m measures the departure 
from the maximally symmetric anti-de Sitter spacetime (see below) and therefore is a 
measure of the inhomogeneity of the spacetime. 

The tetrad components of the fluid velocity one-form of this solution are 



u, 



-1,0,0,0) 



so, the metric is given in comoving coordinates. The non-vanishing Weyl spinor com- 
ponents \& are 



x e 



mS 2 



sinh 



2m+l, 



mS 2 



x)e 



sinh 2m (a;) 



so that the Petrov type is I except when the constant m vanishes and the metric is 
conformally flat. In this case, the density and pressure are p = —3b 2 / 4 , p = 36 2 /4 and 
therefore the spacetime is anti-de Sitter. 
The expansion of the fluid is 



9 = Ssinh 1+m (x) (2-3z) 
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The shear tensor of the solution is diagonal and its components are 

1 c 2S* 2 

an = --V22 = a 33 = - smh 1+m (x), a^a^ = — smh 2+2m (x). 

The acceleration of the fluid is 

a = — cotanh(:r) (m + z) dx. 

The density and the pressure are 

p = S 2 sinh 2m (x) (m - 3z 2 ) , p = S 2 sinh 2m (x) (2m + 3z) z. 

From these expressions it is clear that the perfect fluid does not satisfy an equation 
of state p = p(p) for any value of m. The solution with S being a hyperbolic sine 
has always a region near u — where the energy-density is negative. Thus, we will 
concentrate on the solution with S being a hyperbolic cosine. For this solution the 
energy-density is positive everywhere provided the constant m satisfies 

4 

m>-. 

The range of variation of the coordinates for this solution is 

— oo < t < oo, < x, —oo < y < oo, — oo < z < oo. 

When u — > — oo the function S tends to infinity in an exponential way. So, the length of 
the curve x = const., y = const, and z = const, from any finite value of u to u — — oo is 
finite. Furthermore, the energy-density diverges here and therefore we have a big-bang- 
like singularity at t — — oo. On the other hand, we have that u = +oo is located at an 
infinite distance and the same happens with x — 0, which is also located at infinity. The 
expansion of the fluid is positive everywhere and tends to infinity when we approach 
the big bang and to zero in the infinite future. 

Let us now study in which region the pressure is positive. A trivial calculation shows 
that the pressure is positive in either of the following regions 

. / VI + m 2 \ m-1 ( y/1 + m 2 \ 3 + m 

tanh u > : or tanh u < . 

V 2 j ~ VI + m 2 V 2 / 3^1+^2 

Furthermore, in the range for m we are interested in (m > 4/3), we have 

m-1 , 3 + m 

, < 1 and . > -1, 

yl + m 2 3\/l + ^ 

and, therefore, there always exist two values of u, say m and ui, such that the pressure 
is positive for 

— oo < u < u < and < U\ < u < +oo, 
i.e. it is positive except for a finite interval around u — 0. 
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6.3 Analysis of the solution C 



After rescaling {t, x, y, z] and a coordinate change e" 
takes the form 



x we find that the metric 



ds 2 



-dt 2 + 



dx 2 



S 2 (x,t) 

where the function S(x, t) is 



+ x h{h+1) cosh 1 "^)*/ 2 + x b{b ~ l) cosh 1+b (t)dz 2 



(36) 



i+i>- 



S(x,t) = x 2 + s | sinht 



l+b z 



Sq constant. 



As a switch b <-> — 6 is equivalent to y <-> z, we consider only b > 0. It is convenient to 
define the two functions S x (x) and S t (t) as 



= x 2 , St = so | sinht 
The conformal Killing vector of this solution is 



(S = S x + S t ). 



k = X d x -h(l+b) yd y +h(l-b) Zd z . 



The non-vanishing Weyl spinor components \l/ are 



S 2 



4 cosh t 



b 1-b 2 e 



S 2 



12 cosh^ t 



S 2 



4 cosh t 



-b(l-b 2 ) e*. 



Consequently, the Petrov type of this metric (when 6 7^ 1 and 6 7^ 0) is I everywhere 
except at the spacelike hypersurface cosht =| b^ 1 \ where it degenerates to type D. In 
the particular case 6 = 1, the metric is conformally flat and the density and pressure of 
the perfect fluid are 



p = 3s 



o- 



P = - 3s o> 



so that the metric is the de Sitter spacetime. When b = the Petrov type is D 
everywhere. 

Let us now analyze the kinematical quantities and the matter content for the metric 
(|36l). The tetrad components of the fluid velocity one-form of this solution is 



u, 



■1,0,0,0) 



and the metric is written in comoving coordinates. 
The fluid expansion is 



9 



2 cosh t sinh t 



2S sinh 2 1 - 3S t (l + b 2 ) cosh 2 1 
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The non-vanishing components of the shear tensor and the shear scalar are 

S S 

an = tanht, 022 = — (1 — 3b) tanht, 

3 6 

5 S 2 

a 33 = - (1 + 3b) tanht, a a0 a a/3 = — tanh 2 1 (l + 3b 2 ) . 

6 6 v ' 

The only non-zero tetrad component of the acceleration is 

ai = — S x . 

The expressions for the energy-density and pressure read 

(b 2 - 1) S 2 sinh 2 1 + 3 (b 2 + if S 2 cosh 2 1 



P 



4 sinh 2 1 cosh 2 1 



{b 2 - 1) S 2 sinh 2 1 + (b 2 + 1) S t cosh 2 1 \ 2 (b 2 - 1) S x - (b 2 + 5) S t 
V = 



4 sinh 2 1 cosh 2 1 

The perfect fluid does not satisfy an equation of state p = p(p) unless 6=1 (which 
corresponds to the de Sitter spacetime, as discussed above) or when sq = 0. In the latter 
case the perfect fluid is a stiff fluid (p = p) and the metric coefficients are separable 
functions in comoving coordinates. Therefore, this particular case belongs to the class 
of metrics studied in [|IJ where all the G2 on S2 diagonal stiff fluid metrics with separable 
metric coefficients were considered. 

In order to find the ranges of variation of the coordinates {£, x, y, z} which define the 
spacetime, we must distinguish a number of different cases depending on the signs of 
So and of (b 2 — 1). The reason is that if b 2 — 1 > then the spacetime has no curvature 
singularity at t — and it can be continued across t = to negative values of t. On the 
other hand, when < b 2 < 1 then the energy-density blows up at t = and we have a 
big bang singularity there. Thus, the two possible variation ranges for the coordinate t 
are 

-00 < t < +00, if b 2 > 1 

< t < +00 if < b 2 < 1. 

Regarding the parameter s > we must impose that the function S be non-zero (it can be 
either positive or negative as only S 2 appears in the metric). When sq > 0, x can take 
arbitrarily large values, but when sq < we have that two different variation ranges 
are possible (they correspond respectively to the region where S > and to the region 
where S < 0) . Thus, the different possibilities for the variation range for the coordinate 
x are 

< x < 00, if s > 

2 

< x < (-So) 1 ^ I sinht I if s < 

2 

(—So) 1+62 I sinht |< x < +00 if s < 0. 
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These different variation ranges define completely different spacetimes, they are not 
different regions of a single spacetime. Regarding the coordinates y and z, their variation 
range is, in all cases, 

— oo < y < +00, —00 < z < +00. 

We can now specify which of these different possibilities satisfy p > 0, p — p > and 
p+p > in the whole spacetime. A simple analysis shows that this can be accomplished 
only in one case, namely, when 

b > 1, s < 0, 



2 

< t < +00, < x < (-s )i+^ I sinhi |, 

—00 < y < +00, —00 < z < +00 

(we have to restrict t to t > because of the variation range of x). There is no curvature 
singularity anywhere in the spacetime and so we come to the question of whether it can 
be extended to a larger spacetime or not. The vanishing of S at the boundary 

2 

x = (—So) 1+62 I sinht | 

suggests that the boundary is located at an infinite distance and therefore the spacetime 
can not be extended across this boundary (a formal proof of completeness is more 
complicated, however). Regarding the points with x = the Riemann tensor is regular 
there. The distance to these points along a curve t = const., y = const, and z = const, 
is divergent, which again suggests that they are located at infinity. This does not 
prove, however, that the spacetime is complete although it is an indication for it. Thus, 



this spacetime has no singularities and satisfies the dominant energy conditions |T2 
everywhere. Unfortunately, the strong energy condition p + 3p > is violated in some 
regions of the spacetime. 

6.4 Analysis of the solution D 

The metric is 



-(l+a)r 

ds 2 



e 2ar 

dt 2 + q(r)dr 2 + e~ 2t dy 2 + —-dz 2 

q(r) 



S 2 (r,t) 

where a is an arbitrary parameter and q(r) is given by 



(37) 



aAe 2ar 

q(r) = — - — , A non-vamshmg constant (3c 

1 + Ae 2ar 
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(From the positivity of q follows that A must be negative when a < 0). The subcase 
a = is included by setting A = — 1 and performing the limit a — > which gives 
q = l/(2r). The range of variation of the coordinates y and z are 

— oo < y < oo, — oo < z < oo, 

while the range of variation for the coordinate r depends on the sign of A 

A > -oo < r < +oo, A < r > 0. 

The function S in the metric depends on the variable v = t + r and satisfies the 
differential equation 

(a + 1) S >v (3 + a) S, v 2 S :V S :VV (a + 1) 

H " 7^- H 7^ 1 77- — U. [di)) 



4 S 4 S 2 S 2 I S 

It can be simplified by defining a new function w(v) as 

ciS _ (a + 1) S 



dv 2(w- 1) 



a ^ -1, (40) 



(the case a = — 1 has to be considered separately but we will not study it here in any 
detail as it has negative pressure everywhere). With this equation ( p9|) becomes 



dw w(a + w) . 

T = .I ' ( 41 ) 
at; ty + 1 

which is trivial to integrate (we will not need the explicit solution, though). Further- 
more, S can be integrated explicitly in terms of w from (fH]) as 



S = b- 



in i 1 ~ a 

W — 1 | | U> + a | 2a 

I I 1 + a 

tU 2a 



where 6 is an arbitrary positive constant which gives a global scale to the spacetime. 
This expression for S is also valid when a = by performing the limit a — > 0. The 
result is 



c J ^ - 1 i 
o = o — j j — e 2 » 



The conformal Killing vector of the metric ( [37] ) reads 

k = d t + yd y . 

The Petrov type is D everywhere with the non-vanishing Weyl spinor components \1/ 
given by 
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In the particular case a = 1 the metric is conformally flat and the perfect fluid has 
p = —3b 2 /A and p = 3b 2 /A. Thus, the spacetime in this limiting case is de Sitter 
(A < 0), anti-de Sitter (A > 0) or Minkowski (A = oo, which gives q — 1). 

Returning to the general case a ^ 1, the two repeated null principal directions of 
the Weyl tensor are 



I = 



1 e 



(o+j> 



(dt + e^dy) , k 



1 e 



(o+iv 



The tetrad components of the fluid velocity one-form are 

Ua = , 1 , (--s/g, «>, 0, 0) 



(dt - e _ *dy) . 



and therefore the fluid velocity vector does not lie in the two-plane generated by the 
two repeated null principal directions at each point of the spacetime. 

The condition for the existence of a timelike u is q > w 2 (in the region of the 
spacetime where q = w 2 the energy-momentum tensor represents a null fluid and in the 
region < q < w 2 the matter content is a tachyon fluid, q < is forbidden in order to 
preserve signature). 

The expansion of the fluid is 



e 



iq — w z 



q 2(q-w 2 ) 2 (w - I) 



- 3a(q + w 2 )(q - w 2 ) + q 2 (Aw 2 - 6w - 1)- 
- 2qw 2 {w - l)(w - 2) + 3w 4 



the non-vanishing components of the shear tensor and the shear scalar are 

,2 „„ „ „„2 

O"01 



w 

°oo — — 

9 



w q — w 

pCTllj CT33 + CT22H O"ll=0, 



022 



5 2 /g — w 2 (q + w )(w + 1) 



-e 2 g 



5 + 

r 2 ^ 



(q — w 



2\3 



Iq - w 2 (w - 2q)(w + 1) 



(q — w 2 y 



,2\2 



qw 



The non-zero tetrad components of the acceleration and the acceleration scalar of the 
fluid are 



w 



(l + a)r 

e 2 



Sw< 



a = 



-_d\ 



y/g (g — u> 2 ) 2 (u> 2 — 1) 



2a(g - w 2 ) - q(w 3 + w 2 - 3w - 1) - 2w 3 
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The energy-density and the pressure read 

e (l+a)r S 2 



(q-w 2 



P 
P 



A(w - l) 2 q 

e (l+a)r S 2 



[W 



Ifq 



4q(w - 1) (aw + 2w + 2a + 1) + 3(a + l) 2 (q - w 2 ] 

(a + l)(q — w 2 ) (5aw + w + a + 5) 



q(l - w) (aw + 2w + 2a + 1) - 



4(w + 1) 



It can be easily seen that the perfect fluid does not satisfy an equation of state p = p(p) 
unless a = 1 (i.e. the de Sitter, anti-de Sitter and Minkowski limits of the family, see 
above). The expression for p + p is simple and reads 

p + v S 2 a+a)r (a-w 2 )(a 2 -l) 
2 (1 — w z )q 

and, therefore, in the region where the fluid velocity is timelike we have that 

,2 1 



p + p > 



a 



1-w 2 



> 0. 



In order to study the range of variation of the function w it is convenient to consider 
the metric (pTFD using the coordinates {w,r, y,z}. We will only need the coefficient in 
dw 2 which reads 



-(l+a)r_ 



(w + 1) 



w 



b 2 (w 



w + a 



Tn+rdw 2 . 



Analyzing this coefficient together with the expression for p+p we find that the function 
w is allowed to vary in five different ranges. The extrema of the intervals are 



-oo, 



w 



-1, 



w 



0. 



w 



w = a, 



w 



+oo 



where a may take an arbitrary value. Each one of the five possibilities gives a different 
behaviour for the metric and the matter content. Let us restrict the possibilities by 
imposing p + p > and p > everywhere in the region where the matter content is a 
perfect fluid (that is, in the region q > w 2 ). It turns out that there exist four different 
cases in which these two conditions are fulfilled. The first one is 



a < 



-1 < w < 0. 



In this case both hypersurfaces w — — 1 and w = are curvature singularities of the 
spacetime and they are located at a finite distance. The second possibility is 



a < 



< w < 1. 
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Here, w = is a singularity of the spacetime located at a finite distance and w — 1 is 
at an infinite distance in the future. The third one is 

-f < a < 1, w < -1, 

now the spacetime has a future singularity at w = — 1 at a finite distance and w = — oo 
is at the infinite past. Finally, 

-1 < a < 1, iu > 1, 

which has a curvature singularity at w = 1 at a finite time in the past and extends 
infinitely into the future. 
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7 Appendix 

The aim of this appendix is the analysis of the equation 

+ §S 2 (x) + ^s(x) + f§S 4 (a;) + |e 6 (x) = 0, (42) 

which involves two independent variables, x and t + bK(x). Let us first assume that 
£ 4 = 0. Then the equation reads 

Ei(aO + §S 2 (a:) + ^s(x) + |s 5 (x) = 0. 

If, in addition, S 5 = then this equation is a polynomial of second degree in H/H 
which can be solved (unless the polynomial is identically zero) giving H/H (which is a 
function only of t + bK(x)) in terms of functions of only x. This implies that H/H must 
be constant, thus producing homothetic solutions, against our assumptions. Thus, the 
only possibility is 



S x = S 2 = S 3 = £ 4 = £ 5 = 0, H arbitrary 



which we label as Case 1. 
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When £ 4 = and £5 7^ 0, we can divide equation ([42|) by £ 5 and apply the 
differential operator D = d x — bK'dt which kills all functions depending only on the 
variable t + bK(x). The result is 

Now, the same considerations made before imply that each coefficient in this polynomial 
in H/H must be zero. Thus, we necessarily have 

£ 4 = 0, Si = miE 5 , £ 2 = ^2^5, S 3 = m 3 £ 5 , £ 5 7^ 0, 

ii" ^2 ^ 

mi + m 2 — + ^3-^2" + jj = 0, 

where m 1; m 2 and m 3 are arbitrary constants. We label this possibility as Case 2. 
We have finished the analysis when £ 4 = 0, so let us assume from now on that £4 7^ 
0. Dividing ( f42|) by £4 and applying the operator D on the resulting equation, we 
immediately find 

EiV /£ 2 \'# /£ 3 \'# 2 /£ 5 \' H , . 

1 + 77+ 7H+ ^ 77 = 0- (43) 



.£4/ VS4/ VS4/ i/ 2 VS4/ H 

When (E5/E4)' e0 we necessarily have 

Si = fc]£4, £ 2 = /C2S4, S3 = /C3S4, £5 = /C5S4, £4 7^ 

, , H , H 2 H H H 

h + k 2 h k 3 — H h h— = 0, 

1 2 H 3 H 2 HH °H 

(ki, k 2 , k 3 and k$ arbitrary constants), which is the possibility labeled as Case 3. It only 
remains the situation when (E5/E4)' 7^ 0. We can divide equation ([43]) by (E5/E4)' 
and apply the differential operator D on the resulting equation. We find a polynomial 
involving only H/H which, as usual, implies that each coefficient must vanish. The 
result is 

Ei = fciE 5 + /5iE 4 , £ 2 = a 2 £ 5 + /3 2 £ 4 , £ 3 = a 3 £ 5 + /c 3 £ 4 , 

where k\, k 3 , a 2 , a 3 , j3\ and [3 2 are constants. Equation fl4~3"| ) is 

H H H 
h + a 2 — + a 3 — + — = 0. 
ti ti ti 



It still remains to impose the original equation ( f42|) which, using all this information, 
reads 

H H 2 H 3 

Pi + (P2 - ki) — + (k 3 -a 2 ) — - a 3 — = 0. 



35 



Each coefficient in this polynomial must vanish and we finally find 

Si = &4E5, £ 2 = h^A + faT, 5 , S 3 = /c 3 S 4 , S 5 7^ 0, — ^ const. 

^5 

, H H 

which is the last possibility, labeled as Case 4. 
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